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ABSTRACT 

In the last decade, a number of experiments dealt with the problem of measuring the arrival di- 
rection distribution of cosmic rays, looking for information on the propagation mechanisms and the 
identification of their sources. Any deviation from the isotropy may be regarded to as a signature of 
unforeseen or unknown phenomena, mostly if well localized in the sky and occurring at low rigidity. 
It induced experimenters to search for excesses down to angular scale as narrow as 10°, disclosing 
the issue of properly filtering contributions from wider structures. A solution commonly envisaged in 
these years based on time-average methods to determine the reference value of cosmic ray flux. Such 
techniques are nearly insensitive to signals wider than the time-window in use, thus allowing to focus 
the analysis on medium- and small-scale signals. Nonetheless, often the signal cannot be excluded in 
the calculation of the reference value, what induce systematic errors. The use of time-average methods 
recently brought to important discoveries about the medium-scale cosmic ray anisotropy, present both 
in the northern and southern hemisphere. It is known that the excess (or the deficit) is observed as 
less intense than in reality and that fake deficit zones are rendered around true excesses, because of 
the absolute lack of knowledge a-priori of which signal is true and which is not. This work is an 
attempt to critically review the use of time average-based methods for observing extended features in 
the cosmic-ray arrival distribution pattern. 



INTRODUCTION 

Large field of view (fov) experiments operated for 
cosmic-ray (cr) physics collect huge amount of high- 
quality data, making possible to study the CR arrival dis- 
tribution with remarkable detail. Either satellite-borne 
or ground-based detectors are considered, many collab- 
orations coped with the measurement of the CR inten- 
sity all over the portion of the sky they observed. They 
all looked for deviations from the isotropic distribution, 
as any signature of anisotropy provides essential infor- 
mation on CRs and the medium that they propagate 
through. 

Apart from the search for gamma-ray emission from 
point-like (or quasi point-like) sources, either in the MeV 
energy range on-board satellites and in the TeV region 
with ground-based telescopes, directional data are ana- 
lyzed to map the CR gradient all over the sky at every 
angular scale. 

Signal as deviation from the isotropy 

Any motion of the laboratory system with respect to 
the CR plasma turns to a dipolar signature with a 
maximum in the direction of the motion. This is 
true for any "still system" we want to consider: it 
might be well the Solar System (i.e. the motion 
of the Earth around the Sun is factorized) or the 
Galaxy itself (the motion of the Solar System around 
the Galaxy center is factorized). Such a process of 
dipole generati on is commonly referred to a s "Compton- 
Getting" effect ()Compton and Getting 1935) and was ob- 
served by a number of experiments (jAglietta et al.lll996t 
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lAmenomori et alJ2004l2007HAbbasi et al.l2011f >. In this 
last the amplitude and the phase of the signal 

are analytically predictable, the observation is commonly 
considered as the starting point of any anisotropy anal- 
ysis, as it demonstrates the reliability of the detector 
and the analysis methods to be fine-tuned. Concerning 
the "galactic" Compton-Getting effect, the importance 
of this measurement lays in determining whether the CR 
plasm a is co-moving or not wit h the "still system" under 
study (|Amenomori et al1l2006l ). 

Moving to narrower scales, it is known that a CR "pure" 
anisotropy, i.e. not due to expected Compton-Getting ef- 
fects, exists down to angular scales as wide as ~ 60°. It 
was observed by ground-based detectors ever since 1930s 
and most recent exper iments represented it in 2D sky- 
maps (see lluppal 120121 and references therein) . Such a 
"large-scale" anisotropy (lsa) is of fundamental impor- 
tance, commonly interpreted as a signature of the prop- 
agation of CRs in the local medium. 

Being charged particles, CRs have trajectories deflected 
by magnetic fields, so that their rigidity sets up a "mag- 
netic horizon" , i.e. a distance below which the ob- 
served arrival distribution contains information about 
the interaction of CRs with the medium that they prop- 
agate through. The diffusion approximation effectively 
explains the observations beyond this horizon. Thus, 
GeV-TeV CRs are an effective tool to probe magnetic 
fields within the Solar Sy stem (up to the Heliotail) 
(|Desiati and Lazarianl feOia . The multi-TeV region is 
important to study the CR propagation in the Local 
Inter-Stellar Medium (LISM), whereas higher energy CRs 
may reveal important features of the galactic magnetic 
fields. If electrons (e±) are considered, synchrotron en- 
ergy losses should be accounted for in defining the mag- 
netic horizon (remarkably closer than for protons of the 
same rigidity). Apart from that, the line is the same 
and gives the importance of any attempt to measure 
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the anisotropy even in the e± channel. The CR elec- 
tron anisotro py was recently search ed for by the Fermi 
experiment (Ackcrman n et all 1201 0). though with null 
result. 

Since 2009 "medium-scale" anisotropy (MSA) struc- 
tures were observed in the CR distribution down to angu- 
lar scale as wide as 10°, their origin still unexplained (CR 
source region, magnetic structures focusing CRs, etc. . . ) 
(lAmenomori et all 120071: 1 Abdo et al.ll2008at lluppa et all 
20123). They were observed both in the Northern and 
the Southern hemisphere at TeV energy and they do not 
appear correlated to each other. It is quite natural to 
expect that a well defined CR pattern at any angular 
scale will be found as soon as statistics will be enough. 
The narrower the structures the closer their origin, what 
explains the growing interest towards this phenomenon 
(Desiati and Lazarian l2012tlDrurv and Aharoniarj|l2008t 
ISalvati and Saccoll2008t IGiacinti and Siglll2012D . 

Besides these CR signals, diffuse gamma-ray emission is 
often measured by satell ite experiments in the GeV range 
(|Ackermann et alll2012f ) and extensive air-shower (EAS) 
arrays observed the diffuse emission fro m the Galactic 
Plane at TeV energy (jAbdo et al.ll2008bT ). Structures as 
wide as 10° — 20° have to be properly extracted from 
the background (the overwhelming CR contribution or 
the average photon content) and the analysis methods 
used to do that often rest on the same ideas exploited to 
measure the CR anisotropy. 

Detection techniques 

Either the CR anisotropy or the diffuse emission from 
the Galactic plane is considered, the experimental issue 
of properly detecting and estimating the intensity of a 
signal as bright as 10 -4 — 10 -3 with respect to the average 
isotropic flux of CRs has to be dealt with. 

It translates in estimating the exposure of the detector 
with accuracy well below that threshold, to avoid that de- 
tector effects mimic signals due to physics. To keep the 
exposure map under control down to ~ 10 -4 is a chal- 
lenge even for the most stable experiment, as unavoid- 
able changes in the operating conditions occur and not 
always on line corrections can be readily applied. Both 
for satellite experiments and ground-based detectors the 
envisaged solution is to estimate off-line the exposure, re- 
lying on statistical methods applied to the large data-set 
available. 

Good results can be achieved by combining Monte 
Carlo simulations and t he record of the operati ng con- 
ditions (see for instance ()Ackermann et alll2012D ). Oth- 
erwise, if simulations do not reach the needed accuracy, 
like for EAS experiments^, or just to have a simulation- 
independent result, the exposure is estimated from data 
themselves, by exploiting some (assumed) symmetries in 
the data acquisition. 

A number of data-driven methods to estimate the ex- 
posure exist, although all of them are based on geometri- 
cal properties of the detector a cceptance (see for instanc e 
the "equi-zenith" methods (jAmenomori etaD l2005h l 
and/or on the uniformity of the trigger rate within a 
certain period (Alexan dreas et al 1119931 iFlevsher et all 

3 The important effect of temperature and pressure variations on 
the atmospheric depth and, consequently, on the trigger efficiency 
of EAS arrays cannot be taken into account down to 10 — 4 -10~ s 



120041) . As well pointed out by (jFlevsher et all 120041) . 
these symmetries are assumed to be valid in some con- 
ditions and such an assumption is part of the null- 
hypothesis against which signals are tested. 

Among all the techniques that experimenters devel- 
oped to estimate the exposure, this paper is focused on 
those based on the time-average of collected data. What- 
ever its particular implementation is, any time-average 
method (tam) relies on the assumption that the signal 
to be interpreted as background can be filtered out by 
averaging the event rate in a certain time window. The 
time-average can be performed directly, i.e. by integrat- 
ing the event rate within the time-window and then di- 
viding by the window width (direct integration method), 
or by using Monte Carlo techniques. In the latest case, 
each event is associated with a number N of "fake" events 
having all the same experimental characteristics but dif- 
ferent arrival time, sampled according to the measured 
trigger rate. After "time swapping" (or "shuffling" or 
"scrambling"), the over-sampling factor N is accounted 
for and the final result makes this approach equivalent to 
the Direct Integration method. Actually, the only differ- 
ence is that the integration is performed via Monte Carlo 
instead of directly. 

The use of tams is favored by the property well known 
in signal processing for which smoothing out a signal with 
a top-hat kernel as wide as T strongly suppress signals 
narrower than T. As a consequence, the smoothed signal 
will contain only signals wider than T, so that subtract- 
ing it from the actual signal is the same as saving only 
contributions from frequencies higher than 1/T. As it 
will be discussed in a more formal way in the next sec- 
tions, TAMs demonstrated to effectively work for point- 
like and quasi point-like gamma-ray sources, because in 
these cases a suitable time interval around the source 
can be excluded from the integration. The time average 
alongside the excluded region is quite a robust estima- 
tion of the exposure, i.e. - through a simple scaling with 
the average trigger rate - of the CR background. 

The attention on TAMs recently grew up because of 
their application to detect medium scale excesses on 
top of the large scale CR anisotropy. The Milagro col- 
laboration firstly tried to "adapt" the Direct Integra- 
tion meth od for studies on sm all to intermediate scales 
(10°-30°) (|Abdo et all I2008aft . The attempt came from 
the assertion that averaging for a time interval T cor- 
responds to make the analysis insensitive to structures 
wider than 15° T/l hr in Right Ascension. Afterwards, 
other experiments applied tams for anisotropy studies, 
either to estimate the over all exp osure or to focus the 
analysis on a certain angular scale dGuillian et al 1120071: 
l Ackermann et all 120101 : TAbbasi et all 120111 : lluppa et all 
2012a). In some cases, the property of filtering out larger 
structures became the main reason why tams were used, 
although no detailed discussion was ever made on the 
potential biases of these techniques in filtering the large- 
scale structures. 

This paper collects a series of simple calculations and 
observations on the filtering properties of TAMs. As they 
are applied in a variety of experiments having different 
operating modes, sky-coverage and trigger rate stability, 
a general treatment of the matter is impossible. When a 
specific experimental layout had to be accounted for, the 
authors made use of a virtual EAS array similar to the 
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ARGO-YBJ experiment (jBartoli et all 1201 lal ). whereas 
to discuss a likely case of underlying large scale structure 
the model of the large sca le anisotropy of CRs as given 
in (jAmenomori et al.ll2007f) has been used. 

Statistical effects were not considered, i.e. no Pois- 
sonian fluctuations around the average event content of 
each pixel were accounted for. In fact, this contribution is 
to outline some major potential systematic effects, intrin- 
sic to the application of tams, regardless if the number 
of events is sufficient to make them visible or not. 

The paper is organized as follows. In the section [T] 
an introduction on tams as exposure estimation meth- 
ods is given. The section [5] is a brief interlude which 
demonstrates a consequence of data-normalization along 
the right ascension to be considered for all further discus- 
sions, though mostly affecting the £ = 1,2,3 components 
of the signal. In the section |5] the effect of TAMs on the 
signal to be detected are introduced, mostly for what con- 
cerns the reduction of the intensity and the appearance 
of border effects. The section [4] finally provides quanti- 
tative information on the residual contribution from fil- 
tered components and the signal distortion due to the 
method. Some conclusive remarks are given in the last 
section. 

1. EXPOSURE CALCULATION WITH TAMS 

From the experimental viewpoint, the observation of 
excess (or deficit) effects at a level of 1CT 4 is a difficult 
task, because of the intrinsic uncertainty that CR appa- 
ratus have to cope with in estimating the exposure. For 
EAS arrays the atmosphere is part of the detector it- 
self and data must be handled with care to avoid that 
a atmospheric change mimics a signal somewhere in the 
sky. For detectors on board satellite, no atmosphere ef- 
fects are there but trigger rate variations persist related 
to changing conditions along the orbit. 

In general, assuming there is an isotropic charged CR 
flux overwhelming all the other signals, the exposure is 
estimated by assuming it proportional to the integrated 
CR flux. In this way, the exposure estimation problem is 
posed as a CR-counting problem. 

Hereafter, the number of events collected (or com- 
puted) in the solid angle dfi centered around ft = (9, <fr) 
in the local frame, in the time interval [t, t + dt) will be 
written as: 

dN(Q,t) ( d 2 iV(fi,t)' 



dt 

to lighten the notation. 



dn dt 



1.1. Point-like and quasi- point-like sources 

For point-like or g-uasi-point-like sources, TAMs are usu- 
ally applied to estimate the exposure (i.e. the expected 
background CR rate) from a certain direction of the sky. 
They are an evolution of the elder "on-off" method and 
rely on the assumption that the CR flux from a given di- 
rection Q in the local reference frame is practically con- 
stant during short time-periods. In other words, the av- 
erage count from = (#,</>) during the the interval T is 
quite a good approximation of the CR number N cr : 



dN b (n,t) dN b (n,t) _ / dN ev (n,t) 



dt 



dt 



dt 



(1) 



where N b is the actual (unknown) background CR num- 
ber, N b is the estimated one, and N ev indicates the num- 
ber of measured events. The average is computed in the 
time interval T and using the kernel function uu, so that: 



dN ev {n,t) 

dt 



Jt-T/2 dT W ( T > 



(2) 



If the source contribution is not excluded, the function 
w(t) in ^ is the trigger rate and accounts for over-all 
variations in the acquisition regime: 



dr w(t) = dr 



dn 



dN ev (n,t) 



FOV 



dt 



(3) 



The integration is carried out numerically, with the Di- 
rect Integration or the Time Swapping method (see the 
Introduction). 

The next sections of this paper will focus on the role of 
dN ev /dr in the estimate ©, as this quantity is the sum 
of different contributions and the problem of a proper 
separation of the signal in the angular domain via tams 
has to be approached by considering the time properties 
of dN ev /dr. However, before that, two other aspects of 
the equations ([I])-© should be made explicit. 

• Time interval. The quality of the approximation 
is related to the difference between N b (fl,t) and 

Nb(£l,t) {!]), i.e. to how representative the time- 
average is of each instant CR-flux. If the time win- 
dow T is chosen too wide, the geometrical distribu- 
tion of the CR arrival directions may significantly 
change, due to atmospheric effects. Some changes 
in the detector operating regime may have the same 
effect of making the dN ev /dr distribution not uni- 
form. 

• Source exclusion. The source contribution should 
be excluded from the time-average. Mathemati- 
cally, the weight Q has to be replaced by w se (r): 



w(t) 



if Q G D src (r) 
w(t) otherwise 



(4) 



where D src (r) indicates a confidence solid angle 
around the source at the time r. 

As far as the time window is concerned, the acquisi- 
tion of EAS arrays is not stable for periods longer than 
2 — 3 hrs, as climatic changes affect either the arrival di- 
rection distribution of CR and the detector response to 
the incoming radiation. There are far minor problems 
for underground experiments or neutrino observatories, 
where even longer times ar e used in the literature (up to 
24 hrs (A bbasi et al 2011)). Nonetheless, time intervals 
as short as 4 hrs or less are used also in some of these 
cases to extract small-scale signals. Satellite-borne de- 
tectors usually are so stable to allow to the experimenter 
to shuffle events within the whole data-set available (up 
to few years), so that t he analysis dows not su ffer the 
pitfalls described below (Acke rmann et ail 1201 Of) . 

About the source exclusion, the solid angle to be ex- 
cluded around the source is related to the detector an- 
gular resolution. A safe choice might be 2 or 3 times the 
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average angular resolution plus the source intrinsic ex- 
tension. If a 2°-wide source is observed with an angular 
resolution of 1°, a safe exclusion region of 6° — 8° around 
it can be set. If the region is populated of other known 
sources, the definition of the exclusion region has to be 
obviously adapted. 

For all experiments surveying the sky, the FOV does 
not coincide with the portion of the celestial sphere to 
be investigated. They exploit the rotation of the labora- 
tory frame with respect to the sidereal frame to get the 
project coverage. In this sense, all time-spans may be 
translated into angular intervals measured in the side- 
real frame. If the laboratory rotates around the Earth 
axis (ground-based experiments), time intervals are R.A. 
intervals. Depending on the rotation of the laboratory 
frame in the sidereal frame, 1 hr may correspond to ~ 15° 
in R.A. for a ground-based detector or ~ 240° in the or- 
bit panle of a low Earth circular orbit satellite. For the 
IceCube detector, at the South Pole, the sky portion ob- 
served is always the same and time-flow simply brings 
a rotation with respect to the celestial coordinates. For 
ground-based detectors 2 — 3 hrs correspond to 30° — 45° 
and enough statistics is left to allow the source exclusion 
(~ 50 — 80% of the events inside the time window T can 
be used). In the literature, typical values are found to 
be T = 2 hrs for the t ime interval and A = 6° for the ex- 
clusion region width (jFlevsher et all 120041 ; iBartoli et al] 
1201 lbf) . 

1.2. Wider structures 

If wider structures are considered, the two conditions of 
the previous section cannot be fulfilled at the same time. 
In fact, the off-source integration interval becomes nar- 
rower than the source extension, thus making the on/off 
source event ratio too high and introducing large fluctu- 
ations in the exposure estimation. 

This is true for a number of structures having physi- 
cal extension wider than few degrees. For instance, when 
experiments like Milagro or ARGO-YBJ measure the dif- 
fuse emission from the Galactic Plane, the source exclu- 
sion region is usually a ±5° galactic latitude belt around 
the plane. Studies of systematics are performed by ex- 
tending the region up to ±10°, obtaining a n on-negligible 
contrib ution to the uncertainty (~ 10% (jAbdo et al.1 
l2008bf l). The ±5° choice gives less fluctuations but prob- 
ably still includes some signal events in the background 
estimation. On the contrary, the ±10° is a safer choice 
for what concerns the source exclusion, at expense of the 
statistictQ. Quoting this effect as a source of systematics 
is still acceptable because the experiments do not have 
the sensitivity to extend the measurement up to 10° — 15° 
from the Galactic Plane. Perhaps next-generation exper- 
iments will have it and it will not be possible to exclude 
the whole region of interest when applying TAMs. 

A similar point holds for the MSA regions, often wider 
than 20°, for which the source exclusion is not applicable. 

In these cases, the signal intensity is reduced by a fac- 
tor p depending on the signal and the background mor- 
phology, as well as on the time-window chosen to apply 
the TAM. For uniform background, uniform source with 

4 It should be noticed that ±10° in Galactic Latitude corre- 
sponds to a varying R.A. interval, as the Galactic plane is not 
oriented along the celestial equator. 



extension Tg in local hour angle and time-window T it 
holds p = 1 - T s /T. 

2. TAMS AND LSA 

Before coping with the filtering properties of TAMs, 
we discuss here the effect of TAMs on the measurement 
of the LSA of cosmic rays. Actually, no modern ex- 
peri ment but IceCube us ed TAMs to estimate the expo- 
sure (|Abbasi et al.l [20TlT ) for all-scale analysis, because 
it would mean to average along 24 hrs and to face all 
the issues of detector stability addressed in the previous 
section. Nonetheless, the result reported here is valid 
also for all the oth er measurements of the C R anisotropy, 
e.g. "equi-z enith" (lAmenomori et~aT1[2005r) or "forward- 
backward" (jAbdo et al. 1120091 ) or else. In fact, a common 
device to bypass the ignorance of the absolute detection 
efficiency as a function of the arrival zenith angle (i.e. 
of the declination), is to set the average flux of cosmic 
rays detected in a certain zenith (declination) belt to a 
certain value, the same for all different belts. In other 
words, deviations from the isotropy are not measured 
with respect to the average over the whole sky observed, 
as to do that the efficiency of the detector as a function of 
the zenith must be properly accounted for. Conversely, 
the reference average is computed along each zenith belt. 

We show here that this solution introduces a degener- 
acy in the measurement of the anisotropy, i.e. m = 
components of the signal are suppressed. 

If the signal is looked at as a distribution f(9, <f) on the 
sphere, the act of normalizing the average content of each 
declination belt to zero can be written as the operator 
S: 

f(9, <j>) ^ f'(0, 0) = f{0, eft - i- JJ dcf>f(9, 4>) 

where the /' distribution is the measured one, which 
differs from the "true" / for the average (f)g — 

1/2 J Q 2n dcf> f(9, 4>). We can consider the spherical har- 
monics expansion of the / distribution: 

oo I 

/(M) = E E 4,*r(M) 

£=0 m=-t 

and considering in a closer detail the effect of the average 
on the signal. In fact: 

r #yf(M)= {^(M) mi 

Using the last result, /' can be rewritten as: 

oo 

/'(M) = /(M)-E a °^>>^ ( 5 ) 

where the degeneracy is made explicit. In fact, all terms 
with m = are suppressed by the experimental technique 
applied, what is more important as the multipole order 
I gets lower. 

If the sky is only partially observed, further effects arise 
due to the non-uniform exposure. In fact, if the number 
of events strongly depends on the declination or other 
preferred directions, significant deviations from isotropy 
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(c) 

Figur e 1. Effect of the R. A. normalization on a dipole signal. 
Figure |(a)| represents the input dipole signal, as intens e as 0.1. 
The dipoE vector points at (0 = 63°, = 243°). Figure [(b)] rep- 
resents the dipole reconstructed with methods normalizing the av- 
erage content in each DEC. belt to zero. The dipole vector points 
at (9 = 90°, </> = 243°) and the intensity is 0.089. Figure [(c)] repre- 
sents the difference between the input map and the reconstructed 
one, which turns out to be a dipole as intense as 0.045, pointing at 
= 0°. Notice that 0.089 2 + 0.045 2 = 0.1 2 . 

might be observed only in certain regions of the field of 
view. 

A representation of the effect just described is given 
for £=1,2 in figures QH 

3. FILTERING PROPERTIES OF TAMS 

As time is a synonym for R.A., tams average signals 
along the R.A. direction, i.e. they enjoy the property of 
filtering out large scale contributions to the signal. 

If we have a data series Xk (k = 1, 2, . . . , A/") and we 



(c) 

Figure 2. Effe ct of the R.A. normalization on a quadrupole sig- 
nal. Figure [(a)] represents the input quadrupole signal. Figure [(b) | 
represents the quadrupole reconstructed with metho ds no rmalizing 
the average content in each DEC. belt to zero. Figure [(c)] represents 
the difference between the input map and the reconstructed one, 
which turns out to be proportional to the Yy(8,(f>) function. 

compute for each point the average: 



n+N/2 
k=n-N/2 



(6) 



k = k±N if k < 1 or k > N) then the difference x— £ will 
maintain intact all structures narrower than N , whereas 
all features much wider than N will be suppressed. In 
Fig. [3] we show the results of a toy numerical estimation 
of the time-average effect on the signal intensity estima- 
tion for different angular scales. The red curve clearly 
shows that the average along AT preserves signals on 
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narrower angular scales and strongly reduces wider con- 
tributions. The Fig. [3] trig gers some other considera- 
tions. Firstly the excess (or the deficit) is observed as 
less intense than it really is. This bias can be avoided 
by excluding the source region, what is impossible for 
structures wider than half the time- window extension. A 
second important issue (related to the first one) is that 
fake deficit zones are rendered around true excesses and, 
vice-versa, fake excesses are seen around true deficits. 
The importance of this problem lays in the absolute lack 
of knowledge a-priori of which signal is true and which is 
not. The problem is present mostly for structures com- 
ing from reducing wider features than for really-narrow 
signals. If the actual signal is well separated in the har- 
monic space from the wider structure to be suppressed, 
it may be possible to observe it also without any fil- 
ter, Le^_jus^by_considering the R.A. distribution like 
m (jAbdo et al.ll2f)08afl . On the other hand, if some hy- 
pothesis can be made (from the literature or from in- 
dependent data about the detector exposure), the effect 
of such underlying larger scale signals can be easily esti- 
mated and quoted as systematic uncertainty of the final 
measurement. 

4. TAMS AND MSA 

The study of the MSA in the arrival distribution of CR 
can be approached with different methods. 

Likely the most orthodox way is to evaluate the expo- 
sure with techniques sensitive to any angular scale and 
then apply the spherical harmonics analysis to filter out 
the signal. The use of the a^ m coefficients prevents any 
contamination from harmonic regions other than those 
selected with the (£, m) numbers and allows to define the 
degree of anisotropy in a (mathematically) robust way. 

Another approach, still starting from an all-scale- 
sensitive estimation of the exposure, is to estimate the 
dipole and quadrupole components of the measured CR 
distribution, to subtract them from the experimental dis- 
tribution and to focus on the residuals at scales less than 
90°. For such narrower signals, the analysis is carried 
out in the real domain. 

These two methods enjoy the uncontested feature of 
properly filtering the signal from scales larger or equal to 
90°, although some problems are there due to the par- 
tial coverage of the sphere by the experimental data. In 
fact, either the ag m expansion and the dipole-quadrupole 
determination are achieved with fit procedures over the 
whole sphere: that part of the sky which is empty of 
data has to be suitably masked and the lack of informa- 
tion unavoidably reflects upon the error associated to the 
intensity and the structure position on the sphere. 

Moreover, both the a,£ m and the dipole-quadrupole 
methods rely on an estimation of the exposure all over 
the angular scales, what implies that they can filter out 
the LSA only if it is properly detected. If the all-scale 
analysis revealed some systematics for the LSA, it would 
be difficult to make sense from an ai m expansion of such 
a signal. 

On the other hand, we already hinted that the appli- 
cation of time-average techniques to get the MSA signal 
would introduce systematics on the flux estimation, as 
well as the bias of filtering only along the R.A. direction. 

Nevertheless, two arguments are in favor of these tech- 
niques: 



1 . they do not inherit systematics from effects present 
below the angular scale they are set to filter out. 
In this sense, no Compton-Getting interference is 
expected, neither influence or artifacts induced by 
large scale atmospheric effects, which instead were 
demonstrated to be relevant for the LSA analysis; 
In fact, systematics introduced in misinterpreting 
the detector performance usually affects the whole 
sky hardly being responsible for localized features; 

2. the amplitude of the systematics described above, 
i.e. the residuals from the LSA structures, can be 
evaluated with Monte Carlo simulations if indepen- 
dent measurements or robust models are available. 

4.1. Residuals of underlying large-scale structures 

Before considering how tams reconstruct extended sig- 
nals like the MSA, i.e. which effects the analysis tech- 
nique has on the intensity and the shape of the signal 
under observation, it is worth to investigate the residual 
effect of the underlying large-scale structures, which are 
strongly suppressed by the time-average. 

Signal gradient along the RA direction 

The result of the tam depends on the signal which it is 
applied to, so that no prediction is possible if the signal 
in toto is not considered. Nonetheless, it is convenient to 
focus on two characteristics of the signal separately, i.e. 
the extension and the gradient along the RA direction. 
It is easy to assess that the gradient of the estimated 
background is proportional to the difference of the signal 
gradient at the boundaries of the time average window: 



dt 



dJV e „(fi,t) 



dt 



k I dN ev {n,r) 



dr 



dN ev (n, T ) 



=t+T/2 



dr 



(7) 



=t-T/2 / 



The constant k depends on the kernel function usec0 In 
the numerical implementation, if a "top-hat" kernel is 
used, it turns out to be k = 1. The equality does not 
hold because of the denominator in the equation [3J in- 
troducing second order corrections in dN ev (fi, t) /dt. The 
equation [7] can be figured out by thinking of a discrete 
implementation of the TAM, where a moving time win- 
dow passes from the i-th time-bin to the (i + l)-th. The 
content of the bin centered at ti + \ + T/2 is included in 
the background estimation in place of the content of the 
bin centered at ti — T/2. 

A simple representation of this effect is given in figure 
21 The top panel represent a large scale excess as intense 
as / = 10~ 3 with respect to the isotropic CR flux, drawn 
according to the equation: 

dN ev (n,t) 



dt 

um t-a + Aa/2\ t&nh (t-a -Aa/2 



I 

2 

The signal center is fixed at ao = 12 hrs, the width 
at Aa = 6 hrs and the signal gradient constant 1/w is 

5 It must be symmetric in time, i.e. ui(r) = w(— r). 
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Figure 3. Simulation of the effect of the time average on the estimate of the intensity of signals of different angular extensions. The 
centr al hi ghest signal peaks around 12 hrs, another signal is visible in the spike around 15 hrs and is represented by the blue curve in the 
zoom [(b)] One more signal, around 20 hrs, is represented also in the zoom [(c)] Black curve: sum of all the signals. Red curve: the signal 
of the black curve as it would appear after the application of a 3 hrs-wide TAM. Blue curve: input signal centered at R.A. = 15 hrs. 

changed from oo to 1/1.5 hrs -1 (black to green curves). 
The bottom panel reports how a TAM (top-hat kernel, 
time-window 3 hrs) would filter each signal of the panel 
above. Residuals of the large scale signal remain, whose 
intensity strongly depends on the signal gradient. For 
a non-physical signal like the black one (w — 0), bor- 
der effects as intense as half the input signal are visible. 
These effects are reduced below 10% the input intensity 
if the signal gradient along R.A. is less than 1 hr . It 
can be noticed that residuals are both of positive and 
negative nature. The maximum intensity of the resid- 
uals and their extension (intended as the width of the 
intervals wherein they are always above or below 5% the 
input signal intensity) are reported in the table [1] It can 
be noticed that the width is never above 3 hrs, i.e. the 
time-window used in the analysis. 



Table 1 

Intensity and extension of residuals reported in the bottom panel 
of figure [4] 

The effect of the LSA on the time-average expo- 
sure estimation 

In the last section it was shown how important is the gra- 
dient of the large scale signal intended to be filtered out 
with TAMs. Results were given for a toy-calculation with 
the purpose of enhancing potential biases of the analysis. 
The aim of this section is to evaluate the effect of the 
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Figure 4. Toy-calculation to highlight the relevance of the gradi- 
ent along the RA direction of the signal under study. Top panel: 
a 6-hrs wide signal as intense as 10 - ^ with respect to the underly- 
ing flat background was simulated to rise up with different slope; 
bottom panel: the signal as reconstructed after the TAM-calculatcd 
background subtraction (time- window: 3 hrs). 

residual contribution of a large scale signal whose nature 
is closer to reality. Once again, the result was achieved 
with numerical calculations, and the algorithm applied 
is described in the following. 

lsa parameterization. — To avoid to introduce any circu- 
lar bias, we used the parameteri zation of the LSA given b y 
the Tibet- AS7 collaboration in (|Amenomori et aLll2010l ). 
The authors model their observation with two structures, 
a Global and a Medium Anisotropy. The former signal 
is what is commonly referred to as LSA, whereas the sec- 
ond one lays on smaller scales and is part of the signal 
that the methods discussed in this paper are tuned for. 
It is worth noticing her e that the best fit of the model 
(jAmenomori et al~l [2010). is given after the normalization 
of data along each declination band (see the section [2j . 

Exposure simulation. — The detector exposure was sim- 
ulated by using the local CR distribution obtained for 
a flat ground-based detector, with a standard atmo- 
sphere absorption model (dN/d9 = Iq exp(— kj cos 9), 
with k = 4.8). After time discretization, the local CR 
distribution was computed for each time bin in the side- 
real day and transposed in equatorial coordinates. The 
exposure map showed the characteristic feature of a max- 
imum at a DEC. few degrees above the experiment lati- 
tude (30° N), fading away at the field of view DEC. lim- 
its . As the calculus was performed with events arriving 
within 9 = 50°, the reference DEC. band for this analysis 
is -20° O 80°. 

Event simulation. — As already said, in order to exclude 
uncertainties due to statistics, the event map was not 
filled by following a Poissonian distribution, but using 
the average value expected from the numerical integra- 
tion of the local distribution function. This is the reason 
why no fluctuations are visible in the R.A. profile. The 
actual background map was obtained with the same solu- 
tion: to avoid any fluctuations, any pixel was filled with 
the product of the exposure times the total number of 
detected events. 

Estimated background. — The estimated background map 
was obtained from the event map, i.e. from the sky pic- 
ture containing the LSA signal. The content of each pixel 



was replaced with the average content of all pixel less dis- 
tant than T/2. The TAM was repeated for all values of T 
from 1 to 24 hrs. Results for 1, 2, 3, 4, 6, 12 and 24 hrs 
will be reported only. 

The Healpix "r ing" pixelization scheme was used 
(IGorski et al.ll27)05r i. 

Anticipating the section l4~2| we notice that the angular 
distance between pixels is the same all over the sphere, 
but the R.A. distance increases when the pixel under 
consideration are close to the poles. Consequently, as 
the interval for the average is set in the R.A. space, 
the computation at low DEC. values is carried out on 
more pixels than at higher dec. The methods becomes 
ineffective when the number of pixels at a certain DEC. 
is such that the average along a certain AT is the pixel 
itself, making the background equal to the signal. The 
effect is small for dec. S < 45° and T > 3 hrs = 60°, 
but is important above 8 = 70°, mostly for T < 2 hrs. 

Figure [5] reports the result of the calculus for different 
dec. bands. The plots were obtained by projecting data 
of the event map (dN ev /dt abbreviated with e in the fig- 
ure) and the estimated background map (dN^/dt — > b) 
in the declination interval indicated, then by calculating 
for each bin the ratio (e — b)/b. In every plot, the black 
curve represents the input LSA signal after the R.A. nor- 
malization. As no detector-induced effects nor statistical 
fluctuations are considered, the gray curve perfectly fits 
with it, representing the signal as it would be observed 
with an all-scale sensitive TAM (T = 24 hrs). The other 
curves represent what remains of the LSA structure when 
the T-wide filter is applied. It can be appreciated how 
the LSA is practically suppressed for T < 4 hrs. The ta- 
ble[5]reports the absolute maximum deviation of the LSA- 
induced signal from the zero-reference- value as a function 
of the dec. band and the time-average window, in units 
of 10 -4 . It is worth noticing that the maximum effect oc- 
curs at the maxima of the LSA, so that for medium scale 
structures observed in other regions than those maxima, 
such an effect is far less than reported in the table [3] 



dec. 




Systematic 




band 




error (x 1C 


- 4 ) 






1 hr 


2 hrs 


3 hrs 


4 hrs 


6 hrs 


-20° 0° 


0.07 


0.29 


0.7 


1.1 


2.4 


0° O 20° 


0.06 


0.26 


0.6 


1.0 


2.1 


20° O 40° 


0.05 


0.20 


0.4 


0.8 


1.6 


40° O 60° 


0.04 


0.14 


0.3 


0.5 


1.1 


60° 80° 


0.016 


0.06 


0.14 


0.25 


0.5 



Table 2 

Systematic error induced by the time-average filtering method 
used in the analysis. The LSA was parameterized according to the 
"Global anisotropy" given in (Amcnomori et al. 20TOT). 

Figure [5] makes a point in showing how weak the ef- 
fect of the LSA is when structures narrower than 45° are 
looked for. If the numbers of the tableware considered, 
together with the typical intensity quoted for the MSA 
emission (4 • 10~ 4 — 10~ 3 relative to CR isotropic flux , 
with T < 3 hrs (lAbdo et al.l l2008at lluppa et al.ll2012at 
lAbbasi et al.ll201lD ). systematic residuals from the LSA 
after the TAM are proved to be less than 20%. 

4.2. Reconstruction of MSA structures 
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Figure 5. LSA simulation for different DEC. bands: |(a)| — 20° 0°; 
curves represent different time-average windows (see the legend for 

As the signal to be resolved is not excluded in the 
background computation, distortions will appear in its 
reconstruction. From previous section it should be clear 
that the effect of the tam depends on the actual shape of 
the event distribution, whose composition is not known 
a-priori, i.e. nobody knows what is signal and what is 
background. After considering how the underlying wider 
structures affect the analysis, we describe here the re- 
lation between the actual and the reconstructed signal. 
The results are obtained by considering a top-hat one- 
dimensional signal to which tam with different time win- 
dows are applied. The background is assumed to be con- 



(b)]o° <-> 20° 

ietails). 



(c)p0° <-> 40°;Kd)U0° <-> 60°;Ke^60° <-> 80°. Different 



stant, what is equivalent to assuming that the larger scale 
structures are well separated in the harmonic space from 
the MSA. In this sense, it should be noticed that the 
calculation is the same of that reported in the study of 
the R.A. gradient, with the important difference that the 
ratio of the signal width to the time window was greater 
than 1 there, but it is smaller here. 

The issue of the actual angular size of the signal and 
that of the reduction of the intensity are addressed. 

Signal extension and declination 

If it is true that time is the same of RA, so that time- 
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average corresponds to RA-average, it does not mean 
that the filtering properties of tams are the same all 
over the sphere. In fact, the angle ip between two events 
having coordinates (a.1,6) and (0:2, S) depends on 5: 

cosi/' = sin 2 S + cos 2 5 cos(ai — a%) (8) 

The equation [8] clearly shows the dependence of ip on S 
(as expected, ip = 0° if 5 = 90°). As a consequence, 
the effect of any filter working in the RA space will be 
different according to the declination band: a top-hat 
filter 45° wide in RA, corresponds to a top-hat filter 31.4° 
wide. A representation of the equation [8] for a\ — cti = 
45° is given in figure [6] 



r 50p 




m 



Figure 6. Angular distance between two points on the sphere 
having RA coordinates shifted 45° from each other, as a function 
of the declination. The vertical lines enclose the declination region 
for which many computations of this paper were made. 

This is the reason why tams tend to return structures 
more and more narrow as declination bands farther from 
(5 = 0° are considered. 

Signal reduction 

We consider a signal with intensity dN s /dt and width w 
(in R.A.), above a flat background with intensity dNb/dt. 
It is analyzed with a tam with time window T. The 
(measured) event content is indicated with dN e /dt. The 
reconstructed signal a' 

, _ dN e /dt - dN b ,/dt 
dNb'/dt 

can be compared with the actual one 

_ dN e /dt - dN b /dt _ dN s /dt 
dN b /dt ~~ dN b /dt 

by studying the ratio a = a' /a. The ratio will depend 
on the ratio p = w/T and on a itself. It is easy to 
demonstrate that for < p < 1 and pa <C 1 (conditions 
fulfilled for the measured MSA intensity) it holds: 

a~(l-p)(l+pa)~l-p (9) 

i.e. the reduction of the signal intensity depends linearly 
on the ratio p. 

The equation [5] makes explicit one of the major uncer- 
tainty induced by tams: since the signal width is some- 
thing to be determined, the choice of the time window to 



be used has to be externally directed (i.e. after experi- 
mental trials or other measurements from literature). It 
is worth recalling that the relation [9] is obtained for a 
particular (and non-physical) case and it gives an upper 
limit of the dependence of the signal reduction on the 
quantities p and a. 

CONCLUSIONS 

This paper is intended to contribute to the study of 
the analysis methods implemented to observe extended 
signals in the cosmic ray arrival direction distribution. 
In particular, it focuses on the time-average based meth- 
ods applied when the source cannot be excluded, and 
it points out how they undoubtedly enjoy properties of 
filtering in the R.A. direction, at expense of important 
spurious effects introduced in the signal reconstruction, 
both for what concerns its intensity and its shape. Ex- 
periments like Milagro, IceCube and ARGO-YBJ made use 
of these methods in the last decade (lAbdo et al.l l2008at 
lluppa et alj|2012at lAbbasi et al J 120111 ). 

On the important point of potential residual effects 
coming from larger signals, supposed to be filtered out, 
it has been shown that they depend on the gradient of 
such signals, rather than on their intensity, what makes 
the bias from the known underlying LSA really small. 

Known distortions of the reconstructed signal were an- 
alyzed, giving numerical information about the intensity 
reduction, the presence of border effects (deficits around 
excesses and vice- versa) and a sort of "reshaping" due to 
the filter acting only along the R.A. direction. 

We conclude that the detection of medium and small 
structures with tams may hardly mimic fake signals 
more intense than 2 10 -4 with respect to the average 
isotropic CR flux. Nonetheless, if more detailed studies 
are attempted, like energy spectrum (i.e. signal inten- 
sity) or morphology, care is needed in considering sys- 
tematic effec ts introduced by tams. In the near future , 
experiments (lAbbasi et al.ll2011trAbevsekara et alj|2012t 
iCao et al.ll2010f ) will have the sensitivity to go below the 
10~ 4 level, what will make the choice of tams less and 
less effective. 

In this sense, medium and small-scale CR anisotropy 
are best searched for if a full-scale analysis is applied, 
with standard spher ical harmonics or wavelet techniques 
(|Iuppa et a l. 2012b), allowing the experimenter to focus 
on well-defined regions in the harmonic space. The pos- 
sibility of this approach is related to the capability of 
controlling the detector exposure down to the level of 
the signal to be observed, i.e. accounting for detector 
and environment effects as precisely as 10 -4 — 10~ 3 . Due 
to this reason, time-average methods, if carefully imple- 
mented, can still be considered a good compromise for 
current experiments, which are controlled to this level 
but do not have yet the sensitivity to detect signals as 
low as 10~ 5 . 
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